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We propose an effective lattice Hamiltonian for monolayer M0S2 in order to describe the low- 
energy band structure and investigate the effect of perpendicular electric and magnetic fields on its 
electronic structure. We, first, derive a tight-binding model based on the hybridization of d orbit als 
of Molybdenum and p orbitals of Sulfur atoms and then, introduce a modified two-band continuum 
model of monolayer M0S2 by exploiting the quasi-degenerate partitioning method. Our theory re- 
veals a difference in electron and hole masses and provides trigonal warping effects. Furthermore, we 
predict a valley degeneracy breaking effect in Landau levels. Besides, we also show that application 
of a gate field perpendicular to the monolayer, breaks the vertical mirror symmetry and induces 
different potentials in the three sublayers, can slightly modify the electronic structure including the 
band gap and effective masses. 

PACS numbers: 73.22.-f, 71.18.+y, 71.70.Di, 73.63.-b 



Introduction — Although studies of two dimensional 
(2D) electronic systems go back to some decades ago, 
it was only in 2004 that the first truly 2D one-atom 
thick material, graphene, was isolated successfully [1]. 
Since then the fundamental interest besides the promis- 
ing applications in nanoelectronic devices, has boosted 
the research about atomically thin 2D materials. It has 
been recently demonstrated that monolayer Molybdenum 
Disulphide (ML-MDS), M0S2, a prototypical transition 
metal dichalcogenide, shows a transition from an indirect 
band gap of 1.3 eV in a bulk structure to a direct band 
gap of 1.8 eV in the monolayer structure [2, 3]. The 
electronic structure of ML-MDS exhibits a valley degree 
of freedom indicating that the valence and conduction 
bands consist of two degenerate valleys (if, K') located 
at the corners of hexagonal Brillouin zone. The lack of 
inversion symmetry of ML-MDS results in a strong spin- 
valley coupling and the valence and conduction bands 
can be described by a minimal effective model Hamilto- 
nian with a strong spin-orbit interaction which splits the 
valence band to spin- up and spin-down subbands [4-6]. 
Due to the peculiar band structure, a variety of nanoelec- 
tronic applications [7] including valley tronics, spintron- 
ics, optoelectronics and room temperature transistor [3] 
are suggested for ML-MDS. Induction of valley polariza- 
tion using optical pumping with circularly polarized light 
is validated by both ab initio calculations and experimen- 
tal observations [5, 6, 8, 9]. Also a combined valley and 
spin Hall physics has been predicted as a result of inti- 
mate coupling of spin and valley degrees of freedom [4]. 

In this letter, we propose an effective model Hamilto- 
nian governing the low-energy band structure of mono- 
layer dichalcogenides and show that its electronic proper- 
ties can be tuned by applying a perpendicular gate volt- 
age (PGV). Although our analysis here is focused on ML- 
MDS, our approach can be easily generalized to other 



dichalcogenides systems MX2. We obtain a seven band 
model (ignoring spin splitting) in which four of them are 
contributed mainly from Sulfur (S) p orbitals and three 
remaining mostly originated from Molybdenum (Mo) d 
hybrids. Our theory describes the conduction and spin- 
split valence bands in accordance with early theoreti- 
cal studies [11], recent density functional theory calcula- 
tions [12] and shows energy corrections to the band struc- 
ture by trigonal warping (TW). The physics of naonrib- 
bon, defects, impurities and so on can be studied by our 
lattice Hamiltonian. Intriguingly, our two-band model 
Hamiltonian incorporates new terms when the system is 
subjected to a perpendicular magnetic field. The latter 
gives a Zeeman-like interaction for valley which breaks 
the valley symmetry in contrast to [10]. Next, we intro- 
duce the effect of PGV which leads to shifts in chemical 
potentials of three sublayers consisting of one Mo and 
two S layers. We show that PGV leads to a splitting 
of two high energy bands originating from p orbitals of 
S atoms. One of our main findings is the possibility of 
tailoring the band gap, effective masses and the valley 
splitting of the valence and conduction bands by varying 
the induced potentials in three sublayers. 

Theory and Model — ML-MDS consists of one layer 
of Mo atoms surrounded by two layers of S atoms in a 
way that each Mo atom is coordinated by six S atoms in 
a trigonal prismatic geometry and each S atom is coor- 
dinated by three Mo atoms. Symmetry space group of 
ML-MDS is D\ h which contains the discrete symmetries 
C3 (trigonal rotation), <j v (reflection by the yz palne), 
(Th (reflection by the xy plane) and any of their prod- 
ucts [4]. In addition to the symmetry of lattice, the local 
atomic orbitals symmetries are essential to be consid- 
ered. Trigonal prismatic symmetry dictates that d and 
p orbitals split into three and two groups, respectively: 
{d z *}, {d x 2- y 2,d xy }, {d xz ,d yz } and {p x ,p y }, {p z }. The 
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reflection symmetry along z direction allows the coupling 
of Mo d xzi d yz orbitals with only p z orbital of S atom, 
whose contribution in the valence band maximum (VBM) 
and conduction band minimum (CBM) located at the 
symmetry points is negligible according to the first prin- 
cipal calculations [12]. Therefore CBM is mainly formed 
from Mo d z i orbitals and VBM is constructed by the Mo 
{d x 2_ y 2,d xy } orbitals with mixing from S {p x ,Py} [8, 12]. 

We thus can construct the tight-binding Hamiltonian 
for ML-MDS by using the symmetry adopted states and 
assuming nearest neighbor hopping terms; 

H TB = Y^KA^iv + 4„blb iv + e b ^b%bl} 

+ J2 ^A^hv + b'J+H.c., (1) 

Here a, b(b') indicate the Mo and S atoms in up (down) 
layer, respectively. The indices /i, v show the orbital de- 
grees of freedom labelled as {1,2,3} = {d z 2, d x i_ y i + 
id xy ,d x 2_ y 2 - id xy } and {l',2'} = {p x + ip y ,p x - ip y } 
for Mo and S atoms, subsequently. Therefore the ma- 
trices e a , e b (e b ), and Uj ifJLU =< a, z, ji\H \b, j, v > are re- 
sponsible for the on-site energies of Mo and S atoms, 
and hopping between different neighboring sites in the 
space of different orbitals, respectively. We do need to 
take into account the overlap integrals, <S, defined similar 
to the hopping terms of the Hamiltonian with elements 
Sij^v =< a,i,fjb\b,j,t/ >. 

Due to the trigonal rotational symmetry of Hamilto- 
nian, the on-site energy matrices take the diagonal form 
e c = + U c for c := {a, 6, &'}, in which shows the 
intrinsic value of the on-site energies for the correspond- 
ing orbital state and U c indicates the potential shift in- 
duced by PGV. Moreover, symmetry properties of the 
lattice lead to only three independent on-site energies 
A 1 = ef, A 2 = e% = eg, and B = e\9 = e\9 . Ac- 
cordingly the symmetries imposed by <j v , cr^^C^ result 
in constraints on the numbers of hopping integrals (three 
parameters) and overlap integrals (three parameters) (see 
the supplementary information for details). For the sake 
of definiteness, we choose tu, t 22l and t 2 i as hopping 
integrals and corresponding forms for overlap integral el- 
ements. A good approximation is provided by Slater- 
Koster method [2] in which all of the hopping and over- 
lap integrals are written as a linear combinations of V v d<j, 
Vpdm S p da, and S pd7V where V pda =< R',p, <t|#|R, d, a > 
and Spdcr =< R/,.p, cr|R, d, a >, for instance. To com- 
plete our effective Hamiltonian, we thus need to add 
spin-orbit interaction (SOI) in the model which causes 
spin- valley coupling in the valence band. The large 
SOI in ML-MDS can be approximately understood by 
intra atomic contribution Hso = £(r)S.L. We only 
consider the most important contribution of Mo atoms 
which gives rise to the spin-orbit coupling term H^q = 
A diag{0, s, — s} in the basis of states {1,2,3} where A 



is the spin-orbit coupling and s = ±. To study the 
band structure properties provided by our tight-binding 
model, we find its fc-space form as ^ k ipk s (H — ES)ipk s 

With ^k s = (&ks,l, a ks,2, &ks,3? frks,l', &ks,2'> &ks,l'' ^ks,2 / ) T 

in which c ks ^ = J2i c is^ exp(zk • R*) (c := {a, 6, b'}) are 
the annihilation operator of electrons with momentum k, 
spin s, and orbital degree of p. The Hamiltonian density 
Hks and overlap <Sks are obtained as, 

(H a H t H t \ (IS S\ 

H=\hI H b , S = 5t i o , (2) 
W t H b J \& 1/ 

with the on-site energy Hamiltonian, H a = U a ts + 
diag(Al, A+, H b = (B + U b )t 2 , H v = (B + U b ')t 2 , 
the hopping matrix, 

/tuf(k,uj) -e- <w tn/(k,-o;)\ 
H t = \t 21 f(k,-uj) W(k,0) , (3) 
\-t 22 /(k,0) -e^t 21 f(k : u) J 

and overlap matrix S is defined similar to H t but t^s 
are replaced by s^s. Here, A± = A 2 ± Xs and /(k, gj) = 
e ik-6i +e i(k-5 2 +u;) +e i(k-5 3 -u;) ig the structure factor with 

uj = 27r/3, in-plane momentum k = (k x , k y ), and Si (i = 
1, 2, 3) being the in-plane components of lattice vectors. 

Generally, our tight-binding model leads to seven 
bands for each spin direction, however, in the absence of 
external bias i.e. U c = the symmetry between top and 
down S sublayers reduces the number of bands to five. 
Two of them correspond to the conduction and valence 
bands, from which we calculate the effective electron and 
hole masses, energy gap, and valence band edge. More- 
over, since CBM mostly comes from d orbitals [12], we 
assume 10% mixing with p orbitals for conduction band. 
This provide us with five equations of seven unknown 
parameters based on the values obtained from ab initio 
calculations and experimental measurements. Further- 
more, it is presumable to consider s^/t^ = O.leV -1 
which reduces the number of unknowns parameters to 
five. We consider the energy gap A = 1.9eV, spin- 
orbit coupling A = 80meV, and effective electron and 
hole masses m e = 0.37mo and rri} h = — 0.44rao (mo is 
the free electron mass) [14], and Evbm — — 5.73eV [15]. 
Eventually, all parameters can be fixed and we then ob- 
tain the on-site energies A\ = — 1.45eV, A 2 = — 5.8eV, 
B = 5.53eV and hopping integrals e^^tn = 0.82eV, 
e -W6 t21 = -l.OeV, and e~ i7 </H 22 = 0.51eV. By having 
these parameters, our tight-bonding theory is completed. 

Fig. 1 shows the band structure of ML-MDS consisting 
of five bands for each spin in the absence of external field. 
Two of them are spin polarized (dotted and dashed lines) 
and others are spin degenerate (solid lines). The inset re- 
veals the spin-valley coupling of the valence band while 
conduction band remains spin and valley degenerate. We 
note that due to the limitations of our model, the high 
energy bands may not be comparable with those of first 
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FIG. 1. (Color online) Band structure of ML-MDS consisting 
of five bands in which two of them are spin split in the valance 
band. The dotted lines refers to one spin and dashed line 
denotes another spin component. Solid lines refer to spin 
degeneracy. The inset shows a closer look to the valence and 
conduction bands. Here, ao = acos(0) where a is the length 
of Mo-S bond and is the angle between the bond and x — y 
plane. In the right panel, side and top views of the lattice 
structure are seen where Mo atom (larger green sphere) is 
surrounded by six S atoms. 
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FIG. 2. (Color online) Contour plot of the conduction (left 
panel) and valence (right panel) bands in momentum space 
for spin-up component together with isoenergy lines for eye 
guide. While the conduction band shows almost isotropic 
dispersion, a strong trigonal warping is occurred in the valence 
band around K points and it is due to the difference of the 
orbital structures of the conduction and valance bands. 



principal calculations in a quantitative manner. We fur- 
ther investigate the band structure close to the valence 
and conduction bands via contour plot which shows the 
isoenergy lines (see Fig. 2). A strong anisotropy of the 
constant energy lines can be seen around K points in the 
valence band, due to the TW, while in the conduction 
band all lines are almost isotropic. 

The low- energy Hamiltonian — Here we present an 
effective low-energy two-band continuum Hamiltonian 
governing the conduction and valence bands around K 
and K' points, by exploiting the Lowdin partitioning 
method [17]. We first change our nonorthogonal ba- 



sis (1^)) to an orthogonal one (\i(j f ) = S 1 / 2 ^)), leads 
to a standard eigenvalue problem H\ij)') = E\%j)') with 
H = S~ 1 ^ 2 HS~ 1 ^ 2 . To employ the partitioning method, 
we expand the new Hamiltonian up to the second order 
in q = k — K around K = (47r/3v // 3«o, 0) point which 
can be written as the sum of q independent (Ho) and q 
dependent (Hi) parts, H = Ho + Hi. Then we rotate 
the orbital basis to \ipbi) = 1, - • , 5) which are the eigen- 
states of Hq corresponding to the eigenvalues E^s. In the 
new basis, the transformed Hamiltonian is UqHUo where 
Uo is the unitary diagonalizing matrix. We define two 
subspaces {\ipbi)i ^02)} corresponding to the conduction 
and valence bands and {\ip' 03 ), 1^04)5 1 ^05)} f° r ^ ne three 
remaining bands. Then we take the block diagonal and 
off-diagonal parts of the Hamiltonian at these subspaces 
as i^diag and V respectively and use the unitary trans- 
formation H' = e~°UlHUoe° such that the lowest order 
in V is eliminated. This results in an effective Hamilto- 
nian H' = i^diag + which is block diagonal in the 
subspaces up to the second order in V. 

The final result for the two-band Hamiltonian describ- 
ing the conduction and valence bands reads, 

H TS = ^-a z + Xrs 1 ° z + t a q • a T 

^21 12 

+ -r^M a + faz) + MoQ • <o- x q • a; (4) 
4m 

for spin s = ± and valley r = ± in which cr T = (ra X: a y ), 
q = (qx,q v ), to = 1.68eV, h = O.leV, a = 0.43, 
and j3 = 2.21. Notice that a = rriQ / m + and j3 
///()//// — ArriQV 2 / (A — A) where m± = m e m h/( m h ± m e ) 
and v = todo/h. The Hamiltonian differs from that 
of introduced by Xiao et al [4] because of the second 
order terms in q. The diagonal q 2 terms, which con- 
tribute in the energy, in the same way as first order off- 
diagonal term, are responsible for the difference between 
electron and hole masses recently reported by using ab 
initio calculations [14]. Moreover, the last term leads to 
anisotropic q 3 corrections to the energy which contributes 
in the TW effect demonstrated in Fig. 2. Importantly, a 
vanishes for the case that m e = — m^, however f3 remains 
a constant. 

In order to study the interplay of spin and valley 
physics, we introduce, by ignoring TW, the effect of time 
reversal symmetry breaking term by applying a perpen- 
dicular magnetic field, leading to the appearance of Lan- 
dau levels (LLs) as follows, 

Et TS = ±y [ — 2 — + ^ c ^ n - t )] + 2( ir } 

+ — + fthj c (an- — ), (5) 

where uj c = eB/2mo and Ib = y/frj (eB) are cyclotron 
frequency and magnetic length, respectively. In contrast 
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to the study of [10] , we see a new valley degeneracy break- 
ing term which is the reminiscence of the Zeeman like 
coupling for valleys. As a result, the conduction band 
LLs are valley polarized and the valence band LLs are 
both valley and spin polarized although we have not yet 
considered the usual Zeeman interaction for spins. In ad- 
dition, the n = LL, E^ TS = [A - hu c r(f3 + a)]/ 2 and 
Eq ts = Xrs — [A + hw c r(P — ot)\/2, depend on the mag- 
netic field strength in an opposite way for two valleys. 
More intriguingly, the splitting of LL in the conduction 
and valence bands 5E + « 5AHuj c and SE~ « 4.6fkj c , 
differ with each other due to difference of m e and — m^. 
Furthermore, we can define three extra splitting terms 
for each LL in the conduction and valence bands: val- 
ley splitting 5E+ = E^ s — E+ s w 2.6fkj c , spin splitting 
5Ef = E^_ — E^ + w 2Huj c , and spin- valley splitting 
5E+_ V = £T _ - £~ + « 3.8ficj c . The first and sec- 
ond splitting, present in the conduction band, originate 
from the Zeeman interactions for valleys and spins, re- 
spectively, but the latest comes from both Zeeman terms 
assisting each other in the valence band. 

Incidentally, it is important to investigate the tunabil- 
ity of electronic structure via external PVG. The ver- 
tical bias breaks the mirror symmetry ah and modify 
the on-site energies of atoms in three sublayers of ML- 
MDS. Accordingly, these changes affect the whole elec- 
tronic structure especially the low-energy characteristics 
such as A, ra e , and effective hopping £ - Interest- 
ingly enough, tuned a and f3 due to PVG, can be used to 
emphasis the valley degeneracy breaking when the sys- 
tem is subjected to a perpendicular magnetic field. We 
assume a single-gated device in which the induced po- 
tentials take the values U b = and U b ' = 2U a . The 
variation of mentioned parameters with the induced po- 
tential U b are shown in Fig. 3, where we only illustrate 
the low-energy band structure for different U b values. 
Using simple electrostatic arguments, the induced poten- 
tials for an applied vertical bias V g can be estimated as 
U b = jeVg with 7 w (d/L)e f /e where e, d and e', L 
indicate the dielectric constants and thicknesses of ML- 
MDS and the substrate, respectively. For typical values 
e « 6 [18], e' « 3.9, d « 0.32nra, and L « 300nra with 
Si02 as a dielectric substrate, we obtain 7 w 0.0007. By 
replacing the substrate with high-e gate dielectrics like 
Hf0 2 with e' « 25 [19], the coefficient 7 « 0.004 in- 
creases and leads to U b « 0.4eV for V g = 100V which is 
already used in on-going experiments. Consequently, the 
PGV effects are enhanced by using a proper substrate 
with large dielectric constant. 

It should be noted that the tight-binding model pre- 
sented here with the parameters fitted empirically takes 
advantage of certain simplifications. First, we have ne- 
glected s, p z orbitals of S atom and 5, d yzi d xz orbitals 
of Mo atom mainly because of the symmetry arguments 
and the range of energies which we are interested in. Con- 
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FIG. 3. (Color online) (Left panel) Spin-up conduction and 
valence energy bands, referred to energy E v = Evbm + U/2, 
for different vertical bias values. The conduction band is 
mostly affected by induced potential U — U h . (Right panel) 
Variation of the gap value versus induced potential which is 
almost linear in the chosen range of U. Insets show the change 
in the effective hopping integral and the effective masses with 
U. Using the form of [/-dependence of the parameters, we 
can also estimate the linear change of 8a « —0.15U/eV and 
6/3 w -lMU/eV. 



sequently, our model will be appropriate mostly for the 
low-energy calculations in the vicinity of the conduction 
and valence bands. Another important simplification is 
the ignoring of 12% ionic contribution in Mo-S bonds [16]. 

Conclusions — We formulated a tight-binding Hamil- 
tonian in order to describe the low-energy band struc- 
ture of monolayer M0S2 which can be useful to study en- 
ergy dispersion and transport phenomena in nanostruc- 
tures M0S2. Our model not only successfully reproduces 
the low-energy dispersion of monolayer M0S2, but also 
predicts the difference in the effective hole and electron 
masses and the trigonal warping effect. In addition, the 
modified two-band model leads to a valley degeneracy 
breaking effect in the Landau levels. Finally, we show by 
applying a perpendicular electric field to the monolayer, 
which induces some excess charges in one Mo and two S 
layers that, the electronic structure especially the band 
gap and effective electron and hole masses can be finely 
tuned. 
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SUPPLEMENTARY MATERIAL 
Seven-bands Hamiltonian 

We start constructing an effective tight-binding model 
for a monolayer M0S2 system, assuming the following 
basis orbitals, 

\l) = d z 2, 

12} = -^(\d X Z-y*) + Z|d* y », |1') = -j=(\ Px ) + i\Py)%6) 

|3> = -j=(\ d x*-y*) - i\ d xy)), |2') = -^(\Px) - APv))- 

The Wannier functions for different lattice sites in a 
crystal are localized and ,here, they can be written as 
|c, = cj |0) where i denotes the site and \i indicates 
atomic orbital, c = a or (6, b') shows the annihilation op- 
erator of three different sublattices of monolayer M0S2, 
consisting of one Mo and two S atoms. Up to the nearest- 
neighbor hopping integral, the tight binding Hamiltonian 
can be written as Eq. (1) in the letter, 



= (c,i,/i|#|c,i,z/), 
Uj^v = (a,i,iJ,\H\b(b')J,i/), 
Sij^v = {a,i,n\b(b'),j, v). 



(?) 



The lattice has two important symmetries C3 = 
e -iu;L z /h anc j a ^ . ^ x ^y^ z ) — >. (—x,y,z)} where the first 
one is the trigonal rotational symmetry operator with 
uj = 27r/3 and L z is the z— component of orbital angular 



momentum and the second symmetry indicates the re- 
flection symmetry with respect to y — z plane. The act 
of these symmetry operators on the basis functions can 
be summarized in the following equations: 

C 3 {|2), |2'» = e-{|2), |2')},C 3 {|3), = e"^{|3>, |1')#) 
«t„{|1), |2), |3)} = {|1>, |3), |2)}, a v {\l'), |2')} = -{|2'), 

and C3II) = |1). It should be noticed that, here, we 
have dropped the spin indices. The symmetry relations 
in Eq. (8), impose some constraints on on-site energies 
and hopping integrals, and thus we have 



( A Y 



e a = I A 2 I , e b = e b ' 



B 
B 



(10) 



A 2/ 

hi -e~ lUJ tn \ / e-tn -e-cn 

ts 1± — I hi ^22 , t$ 2 ± = e~ lUJ t2l t 22 



t s 3 



-t 2 2 ~e ZUJ hi 

e~ iuJ hi -£n N 
e %UJ t 2 i t 22 
—t 22 —t 2 ij 



-t 22 -e~^t 21j 



(11) 



Note that the same relations can be found for s§ i± by 
substituting t^s whit s^s. In the presence of spin- 
orbit interaction of Mo atoms, it is easy to generalize e a 
by replacing (e a ) 22 A 2 + \s and (e a ) 33 A 2 -\s. The 
subindexes of the hopping matrices indicate the nearest- 
neighbor vectors, 



y/S 1 

5i± = ^("2" cos ^' — 2 COS ^' ^ sm ^) ' 

5 2 ± = a(0, cos 0, ± sin 0) , 
1 



(12) 



5 3 ± = a{—— cos 



0, — - cos 9, ± sin0), 



where a = 2A3A and = 40.7 degree [1] are the Mo-S 
bond length and the angle between the bond and Mo's 
plane, respectively. To find above equations, we also 
used the operation of C3 and a v on S{± as CsSi± = 5 2 ±, 
C3$2± = #3±> &v$i± = ^3± and (T v 5 2 ± = 5 2 ±. 

Following a method proposed by Slater and Koster [2] 
(SK), all hopping and overlap integrals can be written 
as linear combinations of V p da and V p d^ and S p d a and 
Spdn- In this method, we define some standard hopping 
and overlap parameters as V^/| m | = (R f , /' ,m\H\R,l,m) 

and S^'lml = V 1 m\R, I, m) where m = 0, 1 stands for 
a, 7r bond and other hopping and overlap integrals can be 
found from SK table [2]. In this way we find 

hi = —j= (n x + in y )[V3n 2 z V pd7T + ^(1 - 3n 2 z )V pd<T ], 
1 y/3 

hi = 2 ( n x - in v)[-2 V pd<r(X ~ nl) + V pd7r (l + n%)l}3) 



1 V3 

t22 = -(n x - in y ) 3 ( — Vpda ~ Vpdn), 
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where n = (n x: n y ,x z ) is a unit vector along pointing 
between the nearest-neighbor lattice points. Once again, 
we can obtain similar relations for the overlaps by substi- 
tuting the hopping matrix elements with those of overlap 
matrix. These relations help us to reduce the number of 
independent hopping parameters from three to two. In 
the absence of external bias, due to the symmetry be- 
tween two Sulfur sublattices, we can simply reduce the 



7x7 Hamiltonian to a 5 x 5 as below 

To find unknown parameters, we should first obtain 
the energy bands around the if -point. After solving a 
generalized eigenvalue problem as — ES\iJj} at K- 

point, we can find the energies as 



E x = 
















E 2 = 


(Ax- 


VB)/2- 


185Re(tns^) - 


VUi ^ 


-B)/2-18^e(t u s* u )] 2 ^ 


(1- 


18| Sll | 2 )(18|tn| 2 -^iS) 










l-18|sn| 2 




1 


E 3 = 


(A- 


+ B)/2 - 


- 185Re(t 2 iS2i) - 


- VU- 


+ B)/2-18$te(t 21 s* 21 )] 2 


+ (1 


-18|s 21 | 2 )(18M 2 -A_B) 










l-18|s 2 i| 2 






E± = 


(Ai- 


VB)/2- 


185Re(tns^) + 


VUi ^ 


- B)/2 - mtefaslJ] 2 -f 


(1- 


18|sii| 2 )(18|tn| 2 -AiB) 










l-18|sn| 2 




5 


E 5 = 


(A- 


+ B)/2- 


- 18Ke(t 2 i*2i) + 


- VU- 


+ B)/2- 18^e(t 2 is* 21 )} 2 


+ (1 


-18|s 21 | 2 )(18|t 2 i| 2 -^-B) 










l-18|s 2 i| 2 







(15) 



These eigenvalues include bonding (with the low-energy) 
and anti-bounding ( with the high-energy) of p — d 
bands. Since the conduction and valence bands are 
mostly formed from d z i and d x i_ y i + id xy orbital of Mo, 
therefore, E\ and E 2 are the valence band maximum and 
the conduction band minimum located at if -point, re- 
spectively. 



scribed in the text, we change nonorthogonal basis to 
an orthogonal one and then rotate them by using a uni- 
tary transformation, [To, which diagonalize iio so that to 
reach the new Hamiltonian, 



H 



U$S- 1/2 HS~ 



1/2 



U , 



(16) 



Two-bands Hamiltonian 



Here, we find the low-energy two-band effective Hamil- 
tonian with Lowdin partitioning method [3]. As de- 



where T~L — + %\ and up to the second order in q 
q x +iqy,Ui is given by 



«i(q) = 



2 ri 2 







-3*22*^0 







3t 2 2*qao 
-§e^ 21 *|g| 2 a 



2„2 



ie">tn\q\ 2 G 
3e iuj t 21 qa 
3t 22 q*a 





-3e^t n qa \ 
-3t 2 2<7*«o 



2 r ,2 








(17) 



where ao = acos(#). In order to find the unitary trans- 
formation matrix, eigenvectors of Ho = So~ 1 ^ 2 HoSo~ 1 ^ 2 

— 1/2 

might be obtained first. Fortunately, <S 1 can be ana- 
lytically calculated at K- point, however, for nonzero val- 
ues of q, an iterative method [4] should be used. There- 
fore, in the vicinity of the if -point, we can treat % asa 
perturbation by assuming small enough q values and the 
transformed Hamiltonian, H, can be expressed as the 



#dia 

V = 



sum of two parts where V is a perturbation term, 

V[2x2] 

M3x3] y 
/ii 2 [2x3] N 
v /ii2 f [3x2] 

Now, we can employ the quasi-degenerate perturbation 
theory that is based on the idea of constructing a uni- 



(18) 
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tary operator e in such a way to drop the first-order 
V in the transformed Hamiltonian, H' = e~° He° = 

#diag + V + [ffdiag, O] + [V, O] + ±[[#diag, O], O] + ■ ■ ■ . 

This imposes the constraint V + [Hdi ag ,0] = which 
leads to the following form for the generator of the trans- 
formation, 

/ r?[2x3]\ 
°~ ^t[3x2] J' (19) 



Insertion of matrix form of operators /in, /ii2, and 77 
results in the two-band form proposed in Eq. (4) in the 
text as 

A A 1 - cr z 
H rs = —a z + Xrs — h t a q • cr r 

ft2| Q |2 

+ — ^-(a + /3^)+tiagq.<(7 x q.< (21) 
for spin s = ± and valley r = ±. 



with the property in which 77/122 — ^11^ = ^12- Then 
if' = i^diag + I [V, O] + • • • is an effective Hamiltonian 
with two decoupled subspaces. Following straightforward 
calculations, the effective Hamiltonian of the low-energy 
bands can be obtained as follows 

/Hi-^r/W + ZW}- (20) 
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